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Abstract 

We consider the inverse problem of determining the isotropic inhomogeneous electromagnetic 
coefficients of the non-stationary Maxwell equations in a bounded domain of M 3 , from a finite num- 
ber of boundary measurements. Our main result is a Holder stability estimate for the inverse problem, 
where the measurements are exerted only in some boundary components. For it, we prove a global 
Carleman estimate for the heterogeneous Maxwell's system with boundary conditions. 
Key words: Inverse problems, Maxwell system, Carleman estimates. 

1 Introduction 

In this paper, we discuss the uniqueness and stability in determining the isotropic electromagnetic coeffi- 
cients of the dynamical Maxwell equations, by boundary measurement of their solution. More precisely, 
given a continuous medium with dielectric permittivity A -1 and magnetic permeability occupying 
an open, bounded and simply connected domain !lcM 3 with C°° boundary F = dQ, and T > 0, we 
consider the following problem for the linear system of Maxwell's equations 

D' - curl (pB) = 0, in Q := x (-T, T), 

B' + curl (AD) = 0, in Q, 

V ' (1.1) 
divD = divB = 0, in Q, 

Dxf = 0, B-u = 0, on S := T x (-T,T), 

where the prime stands for the time derivative. Here the electric induction D and the magnetic field 
B are three-dimensional vector-valued functions of the time t and the space variable x = (xi,X2,xs), 
and v = v(x) denotes the unit outward normal vector to T. Moreover we attach the following initial 
condition to (11.11) : 

B(x,0) = Bo(x), D(x,0) = D (x), i£fi. (1.2) 
Assume that \i and A are scalar functions in C 2 (Q) obeying 

n(x) > fiQ, X(x) > An, x G ri, (1.3) 
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for some Ao > and fiQ > 0- Next, in view of deriving existence and uniqueness results for (11.11 ). 
introduce the following functional space 

#(curl,ft) := {u G L 2 (ft) 3 , curlu G L 2 (ft) 3 }, 

and denote by j T the unique linear continuous application from .ff (curl , ft) into // _1 / 2 (r) 3 , satisfying 
7 r ii = u A v when u G (^"(ft) 3 (see [16][Chap. IX A, Theorem 2]). Then, putting 

# (curl , ft) := {it G F (curl , ft), j T = 0}, 

we see that the operator iA, where 

A® := I ° CUrl ^'' > ] , $ = (D,B) G Dom(,4) := i? (curl;ft) x F(curl;ft), 

y —curl (A.) I 

is selfadjoint in 7~L := L 2 (ft) 3 x L 2 (ft) 3 , endowed with the scalar product 



Further, in light of the last line of (11.11 ). set 

#(div0,ft) := {u G L 2 (ft) 3 , divu = 0} and F (div 0, ft) := {u G F(div 0, ft), j n u = 0}, 

where 7 n is the unique linear continuous mapping from H(div ,ft) := {u G L 2 (ft) 3 , divu G L 2 (ft)} 
onto i^ _1 / 2 (r), such that -f n u = u ■ v when u G Co°(ft) (see |[l3[Chap. IX A, Theorem 1]). Since 
Hq := H(div 0; ft) x Ho(div 0; ft) is a closed subspace of % and that C ker A, the restriction 



is, by Stone's Theorem [ 17 ][Chap. XVII A,§4, Theorem 3], the infinitesimal generator of a unitary group 
of class C° in Hq. Thus, by rewriting (ll.ll )- (ll.2l ) into the equivalent form 



we get that: 

Lemma 1.1 Given (Do, Bo) G V there exists a unique strong solution (D,B) to di.il ) starting from 
(Dq, Bq) within the following class 



:= <AD,D) LW + OuB,B) LW , $ = (D,B) G U, $ = (D,B) G W. 



A)^ = -4^0^ : = $ 6 Dom(A)) = Dom(A) n := V, 




with $ = (D, B) J and <£> = (D , B ) 



T 



(D,B) GC°(M;V)nC 1 (M;^)- 



(1.4) 



Moreover it holds true from lfT6ll TChap. IX A, Remark 1] that 



V = i/ ri o(curl , div 0; ft) x H ni0 ( cu rl , div 0; ft) 



where 



i7* o( cur l 5 div 0; ft) = {u G ii' 1 (ft) 3 , divu = and 7*u = 0}, * = r, n. 
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For further reference we notice from Lemma [TTT1 that the solution (D, B) to dl.lMl.2l l actually satisfies: 

(D, B) G n™ =0 C p {[-T,T];Dom(A™~ p ) provided (D ,B ) € Dom(A^) and A, fi € C m (ft), m > 1. 

(1.5) 

The main purpose of this paper is to study the inverse problem of determining the dielectric permit- 
tivity A -1 = A _1 (x) and the magnetic permeability /i -1 = ijl~ 1 {x) from a finite number of observations 
on the boundary T of the solution (B, D) to dl.lb which corresponds to a realistic physical approach. 
This is an important problem not only in electromagnetics (see RTlO but also in the identification of 
cracks/flaws in conductors (see lfl9l ) or the localization of lightning discharges (see [38 ]). On the other 
hand, we obtain a reconstruction result which involves only a finite number of measurements which is 
not the case in most of the existing results. 

1.1 Inverse problem 

For suitable Bq, Dq, k = 1, 2, we aim to determine A(x), /i(x), x € Q, from the observation of 

B*(a?,t), Dj(x,t), (i,t)eE, k = 1,2, 

where B r = B — (B • v)v (resp. Y) v = (D • v)v) denotes the tangential (resp. normal) component of B 
(resp. D). 

Notice that only a finite number of measurements are needed in the formulation of this inverse prob- 
lem. For an overview of inverse problems for the Maxwell system, see the monograph [39 ] by Romanov 
and Kabanikhin. For actual examples of inverse problems for the dynamical Maxwell system involving 
infinitely many boundary observations (this is the case when the identification of the electromagnetic 
coefficients is made from the Dirichlet-to-Neumann map), we refer to Beleshev and Isakov [3], Caro 
lPT2l . Caro, Ola and Salo lfT3l . Kurylev, Lassas and Somersalo [31], Ola, Paivarinta and Somersalo 11331 
and Salo, Kenig and Uhlmann BUI . It turns out that a small number of uniqueness and stability results 
for the inverse problem of determining the electromagnetic parameters of the Maxwell system with a 
finite number of measurements are available, such as IT331I341 . In both cases, their proof is based on the 
methodology of iPTOl or lF23l . which is by means of a Carleman estimate. 

For the formulation with a finite number of observations, Bukhgeim and Klibanov [ 10 ] proposed a re- 
markable method based on a Carleman estimate and established the uniqueness for similar inverse prob- 
lems for scalar partial differential equations. See also Bellassoued (H, (2), Bellassoued and Yamamoto 
B, HI, A. Benabdallah, M. Cristofol, P. Gaitan and M. Yamamoto Q, Bukhgeim HI, Bukhgeim, 
Cheng, Isakov and Yamamoto 0, Cristofol and Roques Ifl4l . Cristofol and Soccorsi lfT51 . Imanuvilov 
and Yamamoto - (231 , Isakov EH, Khaidarov [26], Klibanov (271, EH, Klibanov andTimonov |[29l . 
Klibanov and Yamamoto QUI . Li and Yamamoto |[33l - |[34l . Yamamoto Il43l . 

A Carleman estimate is an inequality for a solution to a partial differential equation with weighted 
L 2 -norm and is a strong tool also for proving the uniqueness in the Cauchy problem or the unique con- 
tinuation for a partial differential equation with non-analytic coefficients. Moreover Carleman estimates 
have been applied essentially for estimating the energy (e.g., Kazemi and Klibanov [25]). 

As a pioneering work concerning a Carleman estimate, we refer to Carleman's paper [11] which 
proved what is later called a Carleman estimate and applied it for proving the uniqueness in the Cauchy 
problem for a two-dimensional elliptic equation. Since iTTTTl . the theory of Carleman estimates has been 
studied extensively. We refer to a general theory by Hormander (201 in the case where the symbol of a 
partial differential equation is isotropic and functions under consideration have compact supports (that 
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is, they and their derivatives of suitable orders vanish on the boundary of a domain). Later Carleman 
estimates for functions with compact supports have been obtained for partial differential operators with 
anisotropic symbols by Isakov [24]. Carleman estimates for functions without compact supports, see 
Imanuvilov CTTl . Tataru Il42l . As for a direct derivation of pointwise Carleman estimates for hyperbolic 
equations which are applicable to functions without compact supports, see Klibanov and Timonov [29], 
Lavrent'ev, Romanov and Shishat-skii l32l . 

The Carleman estimate for the non-stationary Maxwell's system was obtained for functions with 
compact supports, by Eller, Isakov, Nakamura and Tataru lfl8l . Li and Yamamoto 11331 - IT341 . prove 
a Carleman estimate for two-dimensional Maxwell's equations in isomagnetic anisotropic media for 
functions with compact supports. Lemmas 12.11 and 12.21 are our Carleman estimate for the Maxwell's 
system whose solutions have not necessarily compact supports. 

By the methodology by ifTOll or ll23l with such Carleman estimates, several uniqueness and stabil- 
ity results are available for the inverse problem for the Maxwell's system dl.lb . That is, in Il33l - ll34l 
Li and Yamomoto established the uniqueness in determining three coefficients, using finite number of 
measurements. 

Li and Yamamoto [34], consider nonstationary Maxwell's equations in an anisotropic medium in the 
(xi,X2,X3) -space, where equations of the divergences of electric and magnetic flux densities are also 
unknown. Then they discuss an inverse problem of determining the x 3-independent components of the 
electric current density from observations on the plane X3 = over a time interval and prove conditional 
stability in the inverse problem provided the permittivity and the permeability are independent of X3. 

In ll37l . S.Nicaise and C.Pignotti, consider the Heterogenous Maxwell's system defined in an open 
bounded domain. Under checkable conditions on the coefficients of the principal part they proved a 
Carleman type estimates where some weighted i7 1 -norm of solution is dominated by the I? norm of 
the boundary traces d v U and Ut, modulo an interior lower-order term. Once homogeneous boundary 
conditions are imposed the lower-order term can be absorbed by the standard unique continuation theo- 
rem. Unfortunately, to our knowledge, these results may not be applied directly to the linearized inverse 
problem associated to the original problem. 

Our argument is based on a new Carleman estimate. In comparison with [33 ] and (37), our Carleman 
estimate is advantageous in the following two points: 

• We show a Carleman estimate which holds over the whole domain Q. We need not assume that the 
functions under consideration have compact supports and so ours is different from the Carleman 
estimates presented in [33], and we can establish a Holder estimate. 

• We do not need a priori any unique continuation property and compactness/uniqueness argument 
to absorb the lower-order interior term. In our approach, we establish a Carleman estimates for 
^-solutions of the hyperbolic equation with variable coefficients. This is essential to the proof 
of our main result, because here our problem is involved with a source term and we cannot use the 
standard compactness/ uniqueness argument as in ll37l . 

1.2 Notations and statement of the main result 

In this subsection we introduce some notations used throughout this text and state the main result of this 
article. Pick xq G M 3 \f2, set c(x) = /i(x)A(x) for x € cq = fioXo where (jlq and Aq are the same as 



4 



in (1 1 -3b . and assume that the following condition 
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| V log c(x)\ \x — Xq\ < 1 



p_ 



x £ ft, 



(1.6) 



holds true for some p G (0, Co). This purely technical condition was imposed by the method we use to 
solve the inverse problem under study, which is by means of the Carleman estimate stated in Lemma 
12.21 for any weight function ipQ satisfying the two Assumptions (Al) and (A2). More precisely, in the 
particular case where 



ipo(x) := \x — xq\ , x G f2, 



(1.7) 



then (11.61 ) arises from the classical pseudo-convexity condition expressed by (12.51) . The somehow non- 
natural condition (11.61 ) is thus closely related to the peculiar expression (11.71 ) in the sense that another 
choice of -0o fulfilling (Al) and (A2) may eventually lead to a completely different condition on c(x). 

Next, for Mq > and two given functions /J, A" G C 2 (oj), where ui = VLC\0 for some neighbourhood 
O of T in R 3 , we define the admissible set of unknown coefficients n and A as 

K(M ) = { (/x, A) obeying ([OD aiid ||(/*,A)|| ca(n) < M and (p, A) = (/*«, A*) in w} . (1.8) 

Further, the identification of (A, /i) imposing, as will appear in the sequel, that (B, D) be observed twice, 
we consider two sets of initial data (Dq, B$), k = 1, 2, 



D«(x) = (d«(x),d K 2 (x),4(x) 
and define the 12x6 matrix 

f ei x Bl e 2 x Bi e 3 x Bj 



/C(x) 







B fc (x) = ^(x),6^(x),6 3 fe (x) 

\ 



(1.9) 























ei x Dj e 2 x e 3 x 



ei x Bg e 2 x Bg e 3 x B^ 















V 















x G Q. 



(1.10) 



ei x e 2 x e 3 x y 



We then write (Bf (x, £), Df (a;, i)) the solution to (fTTTb with initial data (Bp, Dp), fc = 1,2, where 
(fii, Aj), i = 1, 2, is substituted for (/i, A). 

Finally, noting = H 3 (-T,T; L 2 (T)) n H 2 (-T, T; H 1 ^)) the Hilbert space equipped with 

the norm 

IMI^E) = ll n ll//3(_ TjT;L 2( r )) + |p||if2(_ T)T;J3 -l( r )-j , U G ^(E), 

we now may state the main result of this paper as follows : 



Theorem 1 LetT > c 1/2 max ien \x - x \ and pick (B§, D§) G (H 2 (nf x ff 2 (0) 3 )nV, k = 1,2, 
in such a way that there exists a 6 x 6 minor m(x) of the matrix K,(x) defined in (11.1 Ob . obeying: 



m(x) 7^ 0, x G 
Further, choose (/ij, Aj) G A £J (Mo), i = 1, 2, so that 

< M, fe = 1,2, 



c 3 (-T,T ; iy 2 '°°(n)) 



(1.11) 



(1.12) 
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for some M > 0. Then there are two constants C > and k G (0, 1), depending on f2, w, T, M and 
Mo, such that we have: 



2 



I Mi - /-^Hh 2 ^) + || Ax - A 2 ||H2(n) 

\fc=i 



B 



B* 



+ 



Df-D 



Notice that the condition dl. lib , which is independent of the choice of the unknown coefficients fj, and 
A, actually relates on the initial functions in (1 1 - 2b - Moreover this condition is stable with respect to 
perturbations in C 2 . Namely, if (Bq,Dq) obeys (11.111) then this is the case for (Bq,Dq) as well, pro- 



vided maxfc = i 2 



(Bg,Dg)-(B fc ,D fc ; 



C(«) 



is sufficiently small. Furthermore there are actual choices 



of (Bq, Dq), k = 1, 2, satisfying (11.11b . This can be seen by taking 



Bjjj(x) = ei, Dj(x)=e 3 , B^(x) = e 2 , Dg(x) = e 2 , x G O \ 



and selecting the 6x6 minor formed by rows 2, 3, 4, 9, 10 and 12. 

Theorem Q] asserts Holder stability in determining the principal part within the class defined by (11.8I ). 
under the assumption dl. 12b . Notice from (11.5b that such a condition is automatically fulfilled for A, /i G 
C 7 {Vl) by chosing the initial data (Bq, Dq), k = 1, 2, in Domf^) (which is a dense in 
The proof of Theorem[TJis based on a Carleman estimate stated in Lemma l2~2l under the conditions (11.3b 
and (11.6b . Notice that (11.6b . which is essential to our argument, is much stronger than the usual uniform 
ellipticity condition. 

The remainder of the paper is organized as follows: a Carleman estimate for the Maxwell system 
(ll.lb is established in Section 2, while Section 3 contains the proof of Theorem [T] 



2 Carleman estimate for Maxwell's system 

As already mentioned, this section is devoted to the derivation of a global Carleman estimate for the 
Maxwell system (ll.lb . 



2.1 The settings 

Let us consider the following second order hyperbolic operator 

Pu = dtu(x,t)-di\(c(x)Vu)+& 1 (x,t;d)u, t G R, 

where M\ is a first order partial operator with L°°(f2 x R) coefficients, and c G C 2 (Q) obeys 

c{x) > Co, x G 12, 

for some positive constant cq. Putting 

a(s,f) = c(x)|£| 2 , i£fi, £GK 3 , 
and recalling the definition of the Poisson bracket of two given symbols p and 



(2.1) 



(2.2) 



(2.3) 



dp dq dp dq / dp dq dp dq 
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we introduce two assumptions. 

Assumption (Al). There exists tpo G C 2 (Q; R+) satisfying 

{a,{aM}(x,t)>0, xcTi, ?el 3 \{0}, (2.4) 

where a is given by (12.31) . 

Since Q is compact and a(x,£) is a homogenous function with respect to £, it is clear that (12.41 ) yields 
the existence of some constant g > such that we have: 

^{a,{a,^ }}(x,0>2gc(x)\C\ 2 , leH, £gR 3 \{0}. (2.5) 

Assumption (A2). The function iftofa) nas no critical points on S7: 

min |V^o(^)| 2 > 0. 
xefi 

Further, g being the same as in (12.51 ). fix 5 > and /3 G (0, g), in such a way that, upon eventually 
enlarging T, we have: 

f3T 2 > maxip {x) + 5. (2.6) 

Hence, picking /3q > and setting 

V>(x,t) =Vo(^)-/3t 2 + /3o, ^gH, t€ [~T,T], (2.7) 



so that 



min ip(x, 0) > /3o, 



we check out from (12.61) that 

if)(x, ±T) </3 -5, x G H. (2.8) 

Notice from (12/71 ) and (f2T8l > that 

maxV(x,t) </3 --, |i| G (T — 2e, T], (2.9) 
sen ^ 

for some constant e G (0, T/2). 

In view of (12.71 ) we may now recall the following global Carleman estimate for second order scalar 
hyperbolic equations, with weight function ip : SI x R — > R defined as 

<p(x,t) = e^ x ^\ x G Q, t G [-T,T], (2.10) 

for some fixed 7 > 0. 

Theorem 2 Assume (A1)-(A2). Then there exist two constants Co > and sq > sucn that for every 
s > sq the following Carleman estimate 

C I e 2stp s (\Vv\ 2 + \d t v\ 2 + s 2 M 2 ) dxdt 

< J e 2sip \Pv(x,t)\ 2 dxdt + j se 2stp [{\Vv\ 2 + \d t v\ 2 ^j + s 2 \v\ 2 ^j dadt, (2.11) 

holds true whenever v G H 1 (Q) verifies d-[v(±T, •) = for j = 0, 1, and the right hand side of d2. 1 1 D 
is finite. Here P is defined by it2.iD - it2.2D and At denotes the volume form ofT. 
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For the proof see Bellassoued and Yamamoto [6], where this result is obtained from a direct computation 
based on integration by parts. 

It is worth mentioning that there are actual examples of functions i/jq fulfilling (A1)-(A2), provided 
the conductivity function c defined in (12.11) verifies (11.61) for some xq G M 3 \f2 and g G (0, cq), where Co 
is the constant defined in (12.2b - Indeed, by putting if}o(x) = \x — xq\ 2 and recalling (I2.3I) . we get through 
an elementary computation that 

~ {a, {a, V>o}} (x, = 2c 2 (x) (l - V - ' ( *~ Xp) ) | £ | 2 + 2c(Vc • • (x - x )), 



so (I1.6I ) immediately yields 

i{a,{a^o}}(x,0>2^c|e| 2 - 
This entails (Al) by (I2.2I) . Moreover (A2) is evidently true as well since Vtpo(x) 7^ for every x £ Q. 

2.2 Decoupling of the system of equations 

Consider now the following Maxwell system 



U' - curl (/^V) = f , in Q, 

V + curl (AiU) = g, in Q, 

div U = div V = 0, in Q, 

U x 1/ = 0, V • 1/ = 0, on S, 



(2.12) 



where the source terms f, g G H l {Q; M 3 ) satisfy the boundary condition 

f (x, t) = g(x, t) = 0, (x, i)£wx (-T, T). (2.13) 
For further reference we recall from (12.91 ) that 

• = P 7(A>-5/2) |^p[r_^ 

(2.14) 



max i6 ^(i,i) < d := e^"" 5 / 2 ), |i| G [T- 2e,T), 
min^gQ ¥>(x,0) > di := e 7ft) , 



and then state the main result of Q2.2 



Lemma 2.1 Assume (A1)-(A2) and let h = (f,g) G ^(QjM 3 ) 2 obey dZTTI) . 77ien we may /wid fwo 
constants C\ > and Si > 0, for which the Carleman estimate 

dj e 2s(p s (j V^W] 2 + s 2 |W| 2 ) dxdt < J e 2s ^ (^\V x , t h\ 2 + \h\ 2 ^dxdt + sh 2d ° s \\W\\ 2 Hl{Q) 

+ J se &v (|VW| 2 + |W / | 2 + s 2 |W| 2 )d<7di, (2.15) 
is true for any W = (U, V) solution to the Maxwell system \2.12\ . whenever s > s\. 
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Proof . The first step of the proof involves bringing (12.121 ) into two independent systems of decoupled 
equations. Namely, by differentiating the first line in (12.121 ) with respect to t, and then substituting 
g — curl (AiU) for V' in the obtained equality, we obtain that 

U" + curl Oicurl (AiU)) = f ' + curl (fug), in Q. 

This entails U" + curl (/xiAicurlU) + curl (/xiVAi x U) = f + curl (/iig), and hence 

U" + /iiAicurl (curlU) + V(^iAi) x curlU + curl (//1VA1 x U) = f + curl (//ig), in Q. 

From this, the well-known identity curl curl U = —A U+ Vdiv U and the third line of (12.121 ) then follows 
that 

U // -/xiAiA(U)+^iU = f + curl(/iig), in Q, (2.16) 

where M\ = M\(x, d) is some first order operator with bounded coefficients in Cl. 
Arguing in a similar way, we find that 

V" - ^iAiA(V) + y x V = g' - curl (Aif), in Q, (2.17) 

for another first order operator S\ = S?\{x, d) with bounded coefficients in Cl. 

Therefore, putting (12.121 ) and (I2.16I )- (I2.17I ) together, we end up getting that any solution W = (U, V) to 
the Maxwell system (12.121) verifies 

f U" - /xiAiA(U) + ^iU = d, in Q f V" - /^iAiA(V) + = G 2 , in Q 

< and < 

[ U x v = 0, curl (AiU) • v = 0, on S, [ V • v = 0, curl (/jiV) x z; = 0, on S, 

(2.18) 

where G\ = f + curl (/iig) and G*2 = g' — curl (Aif). 

Further, consider a cut-off function 77 G C°°(M; [0, 1]) fulfilling 

f 1 if It! <T-2e, 
77 t = <^ (2.19) 
[0 if \t\ >T — e, 

where e is the same as in (12.9b . and set 

U tt = r;U, V s =r/V, #i = r?Gi + 2rf\f + r/'U, K 2 = r)G 2 + 2 V 'V + V "V, 
in such a way that we have 

f Uf - MiAiA(U B ) + ^Uj = ifi, in Q f Vjf - /iiiAiA(Vjj) + ^V tt = # 2 , in Q 

[ Ujj x z/ = 0, curl (AiUjj) • z/ = 0, on £ 1 Vj • v = 0, curl (miVj) xz/ = 0, on S, 

directly from (12- 18b - Then, each of the two above systems being a principally scalar hyperbolic system, 
it follows from the two identities 

U tt (-,±T) = UJ(-,±T) = 0, V S (-,±T)=V' S (-,±T) = 0, (2.20) 

and Theorem that Wj = (Ujj, Vg) obeys the Carleman estimate 

C I e 2sLp s (|V Xit W s | 2 + s 2 |W S | 2 ) cfecft 
< V / e 2sv, |l£j(x,t)| 2 dsdt+ / se 2sv (|VW tt | 2 + |W'J 2 + s 2 |W S | 2 ) Ardi, (2.21) 
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for all s > so- Here we have used Theorem[2]for the diagonal system U'l — jUiAiA(Ujj) and that we can 
absorb the non-decoupled first term ffl\XJ». 

Moreover, as rf and rj" both vanish in (— T + 2e, T — 2e) by (12. 14b . there is a constant C > such that 



[ e 2sv \Kj(x,t)\ 2 dxdt<C f e 2s ^ (\h\ 2 + \V x , t h\ 2 ^dxdt + e 2d ° s 



HHQ) ' 



3=1,2 

according to (I2.20I )- (I2.21I ) and since 



J e 2s<p s (| V^tWI 2 + s 2 |W| 2 ) dxdt<C J e 2s ^s (\V x , t W$\ 2 + s 2 \W$\ 2 ^J dxdt+s 3 e 2dos \\W\\ 2 H 



'Q 

we obtain the result. 



HQ) 
□ 



2.3 Reduction of the boundary terms 

The method used to derive a global Carleman estimate for the solution to (11.11 ) is to replace the local 
boundary problem (11.11 ) in Q, x (— T, T) by an equivalent one stated on the half space x (— T, T). 
This is possible since the boundary F can be represented as the zero level set of some C°° function in 
R 3 . Namely, F being a C°° surface, there exist 6 € C°° (R 3 ) and some neighbourhood V of T in R 3 such 
that T = {x G V, 6>(x) = 0}. We choose V so small that V C O, where O is defined in CT2l write 
2/ = (2/1,2/2,2/3) = (2/', 2/3) tne system of normal geodesic coordinates where y' = (7/1,7/2) are orthogonal 
coordinates in F and y% = 6*(x) is the normal coordinate, and call x = $(?/), where $'(y) > for all 
y G V := <I> _1 (V), the corresponding coordinates mapping. As 



f : =r 1 (r) = {!/€V; 2/3 = o} C 



we may assume that V = X (V) is a cylinder of the form T x (— r, r) with r > 0. 
Further, the Euclidean metric in R 3 inducing the Riemannian metric with diagonal tensor g, 



g{y) = *$'(y)$'(y) = Diag(#i, g 2 , 53), 2/ € V, 

we use the notations of [ 36 1 and note — h= 7P- , — h= jp- , — h= tP- the orthonormal basis associated by q 
to the differential basis of vector fields gj^-, ^ . For any vector field X(x) expressed with respect 
to the Euclidian basis -jP- , tP- , - 
X(y) with respect to the new basis vectors 



as X(x) = Yli=i al ( x )'8§r^ we nave an alternative representation 



1_ d 
9i dyi 



1_ d 
92 dy 2 ' 



1 d 

dy 3 



, given by 



X(y) = Y^a\y)- 
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i=i 



9i dyi 



a(y) = **(i/)a(*(y)), tffo) = *'(y)r 1/2 (l/), 2/ € V. (2.22) 



The divergence (resp. curl) operator of any vector field X(y) = Ya=i ^(^TTjf z|r w ^ res P ec t to the 
local coordinates (yi , y 2 , 2/3) is denoted by div g X (resp. curl g X), and can be brought into the form 



div ff X 



1 3 1 



\/det 5 



i=i 



9j % 



Vdety S J (y) 



3 

£ 

3=1 



1 <9c^' 
/9j 9yj 



+ A(y)-X(y), 



(2.23) 
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where A(y) is some three dimensional vector (resp 

3 



1 dff 



1 <9S-M 1 3 1 <9 



5i % / 



(2.24) 



where M(y) is some matrix function), whereas the outward normal vector field to T at y G T is given by 

1 d 



53 5y 3 



(2.25) 



In light of (I2.22I) - (I2.25I) . we find out by performing the change of variable x = <&(y) in (12.121) . the space 
variable x being restricted to be in Q n V C u n V, that 



(2.26) 



U' - curl a (^iV) = 0, in V x (— T, T), 
Y + curl 9 (AiU) = 0, inVx(-T,T), 
div g U = div g V = 0, mVx(-T,T), 
Ux? = 0, Vu = 0, on f x (— T, T), 

where we have set 

£i(l/) = **(l/Vi(*(l/))*(y), My) = **(y)Ai($(y))*(y), yeVcrV). 

Here we used the identity h($(y),t) = (f($(y),i), g($(y), *)) = for i € (— T,T) and y G V, arising 
from (12- 13b - Further, noting U = (Si, S2, S3) and V = (vi,V2,Vz), the last equation in (12.261 ) reads 
Si = U2 = S3 = on T, so we find that 



U f = 0, Vp = 0, on r, 



(2.27) 



where Uf (resp. Vp) denotes the tangential (resp. normal) component of U (resp. V). From this, (12.231 ) 
and the third line in (12.261 ) then follows that 



1 du a 



whence 



93 dy 3 
VS 3 ^ 



1 9Si 1 duo 
H — — - 

92 dy 2 



9i dyi 



A(y)-V = -A(y)-V, 



ei + -^7^e 2 - A(y) • U e 3 , on r. 



51 



92 dy 2 



(2.28) 



Further, as VSi = -J=^e 3 and VS 2 = -|jfff e 3 on f , according to d2~!27l) . dX28l) then yields 



VU 



£|VS/<C 

3=1 



dy-3 



+ 



du'} 



dy3 



+ 



dyx 



+ 



dy 2 



Ur 



on r, 



(2.29) 



where, for the sake of notational simplicity, we shall use the generic constant C > in the remaining of 
§2.31 On the other hand, since 



curl g \J xV 



1 dui 1 duz 



1 







1 <9Sq 



1 du,2 



93 dy 3 y/g\ dyx J ^/gm %l V V^2 dy 2 y/g~3 dy s J ^9293 892 



1 



8 



+ M(y)XJ 9 , 
(2.30) 
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by d2.24b - d2.25l l. we have 



dy3 



+ 



dy3 



< C 



curl g XJ 



+ 



In view of (12.29b . this entails 





2 






2 


vu 


< 


°( 


curl g \J 


+ 










2 




< 


H 


curl 9 U 


+ 



Ur 



+ 



du 3 



dyi 



+ 



du 3 



dy2 



, on r. 



dyi 
V ? U 



+ 



+ 



dy2 



+ 



, on r. 



As a consequence we have 



whence 



VU 




curl 9 (AiU) 


2 

+ 




2 

+ 





vu 



< c 



V 



+ 



v f u r 



+ 



, onT, 



by the second line of (12.26b . 

Similarly, as div g V = in V from the third line of (12.26b . we get from (12.23b that 



1 dv 3 



This, combined with (12.27b . yields 



1 dvi 1 (9wo , 
5i etyl V^2 <9y2 



|W 3 | 2 < C 



dv-i 



and consequently 



VV 



£|V%| 2 <C 



9yi 



V ? V 7 



+ 



+ 



+ 



+ 



dy3 



+ 



dy 3 



(2.31) 



, on T. (2.32) 



Furthermore, in light of (12.27b and (12.30b where V (resp. dj, j = 1,2,3) is substituted for U (resp. 



Uj, j = 1, 2, 3), we see that 



dvi 
dya 



+ 



dv2 
dya 



is upper bounded, up to some multiplicative constant, by 



curl gV 



+ 

follows that 



V? 



and hence by 



curl 9 ( / uiV) 



+ 



V? 



, on T. From this and the first line of (12.26b then 



dy3 



+ 



dvo 



dy3 



< C 



U' 



+ 



v 9 



so, we end up getting with the aid of (12.32b : 



VV 




U' 


2 

+ 


V 9 % 


2 

+ 





onf, 



onT. 



(2.33) 



Finally, putting (12. 15b . (12.3 lb and (12.33b together, we may state the main result of 
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Lemma 2.2 Assume (A1)-(A2) and put h = (f, g). Then there are two constants C2 > and S2 > 
such that the following Carleman estimate 

C 2 s j e 2sip (s 2 |W| 2 + IV^WI 2 ) dxdt 
< f e 2s * (|h| 2 + \V X ,M 2 ) dxdt + ^(W) + s 3 e 2d ° s ||W||^ 1(Q) , 

where 

^(W) = Jse 2stp (|V r V T | 2 + |V T U„| 2 + |U;| 2 + |V;| 2 + s 2 (|U^| 2 + |V T | 2 )) dadt, (2.34) 
holds true for every solution W = (U, V) to H2.12i , provided s > S2- 

3 Inverse problem 

This section contains the proof of Theorem [T] which is divided into five steps. Firstly, the unknown 
parameters A and \i are brought to the source term of the linearized system associated to (ll.ll ). governing 
the variation induced on the solution to (11.11 ) by perturbating the permittivity by A and the permeability by 
fi. The second step follows the idea of Bukhgeim and Klibanov presented in iTTOl . which is to differentiate 
the linearized system with respect to t in order to move the unknown coefficients in the initial condition. 
The next step is to bound the energy of this system at time t = with the aid of the Carleman inequality 
of Theorem [2] The fourth step involves relating A and fi to the above mentioned estimate through the 
Carleman inequality for stationary (div, curl)-systems, stated in Lemma |3~T1 This is rather technical and 
lengthy so we proceed in a succession of the two Lemmas 13.31 and l3~4l The last step, detailed in £13.31 is 
to derive the desired result from the estimates established in Lemmas I3.3H3.4I 

In the remaining of this text, x$ is a fixed point in M 3 \f2, we choose as in (11.71 ) i/jo(x) := \x — xq\ 2 for 
every x G Q, and (po(x) := ip(x, 0), where ip denotes the function defined by (12.71) and (12- 1 Ob - Moreover, 
for the sake of notational simplicity, we shall use the generic constant C > in the various estimates of 

gnu 

3.1 Linearized inverse problem 

Given (^,,Ai) G A W (M ), i = 1,2, and (Bg,Dg) G H 2 (n) 3 x H 2 (n) 3 , k = 1,2, we consider the 
solution (B^, D^) to the system (11.11 ) where (Aj, //j) is substituted for (A, /x), with initial condition (11.21 ) 
where (Bo, Do) = (B§,Dq). Hence, putting 

M = Ml - M2) A = Ai - A 2 , 

and setting 

f fc = curl GuBl), g fc = -curl(AD^), (3.1) 
we find by a straightforward computation that Ufc = — and Vfc = B^ — B2 satisfy the system 



u fc - curl (MiVfc) = f fc> in Q, 

V' k + curl (AiU fc ) = g k , in Q, 

divUfc = divVfc = 0, in Q, 

Ufc X v = 0, Vfc • v = 0, on S, 



(3.2) 
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with the initial data 

U fc (x,0) = 0, V fc (x,0) = 0. (3.3) 
Further, by using the following notations 

X kjj (x,t) = d{x k {x,t) for X = U,V,f,g and j G N*, (3.4) 

it turns out by differentiating (13.21 ) j-times with respect to t that 



U fc,j ~ curl (^ lV fc-i) = f k,j, in Q, 

Y k j + curl (AiU fcj ) = gk,j, in Q, 

divUjtj = div Vfc j = 0, in Q, 

Ufe j x v = 0, Vfcj • ^ = 0, on S, 



(3.5) 



and, due to (13.1I) - (13.2I ). that \J k ,i and Vfc i satisfy the initial condition: 

U M (ar,0)=carl( M Bg) J V M (x, 0) = -curl (AD*). (3.6) 

As will appear in EB.2H3.3I the main benefit of dealing with (I3.1I )- (I3.6I ) in the analysis of the inverse 
problem of determining A and \i, is the presence of these two unknown coefficients in the initial condition 
(T331). 



3.2 Preliminary estimates 

Let j and k be in {1, 2}. As W k j = (U k j,Y k ,j) is solution to (13.51) . we notice from Lemma I2T21 that 

C 2 s f e 2sLp (s 2 |W fcj | 2 + \V x , t W k>j fjdxdt 

< (|h fc / + |V^h fej f) dxdt + @sA W k,j) + s 3 e 2d0S ||W fcii ||^ 1(Q) := 3fcj(s), (3.7) 

for every s > s 2 , where h k j = (f k j, g k .j) and 9§ a ,tp is given by (12.34b . 
Moreover by the assumption (11.111 ) we can derive that 

2 2 2 | 2 

B§(aO > c *> and ^ M$(x) > c*, x G 0\w (3.8) 
fe=i fe=i 

for some positive constant c* . Indeed, ifBj^x) / OorB^x) 7^ for all a; G then ^ 2 =1 |Bq(x)| 2 > 
in the compact set Now if there exist x\ G 0\w such that Bq(xi) = then by (11.111 ) we have 
Bq(x) ^ for any x G 

Further we recall from [36] the following Carleman estimate for stationary (div, curl)-systems: 

Lemma 3.1 There exist two positive constants S3 and C3 depending only on i^q and £1, such that we 
have 

C 3 s [ e 2sip0 |u| 2 dx < [ e 2sip0 (\curlu\ 2 + \divu\ 2 ) dx, (3.9) 
Jn Jn ^ ' 

for every s > S3 and u G Hq (Q). 
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Prior to the proof of Lemma [331 we establish the following technical result, needed in the derivation 
of (l3~T0b - (l3~TTT) . 

Lemma 3.2 There exists a constant > depending only on T such that we have 
\z(x,Qt)\ 2 dx < 2 ( s / \z(x,t)\ 2 dxdt + s^ 1 / \z'(x, t)\ 2 dxdt 

V Jq Jq 

for alls > s* and z G H l (—T, T; L 2 (Q,)). 

Proof . Let r? G C°°([-T,T]; [0, 1]) fulfills (|2~T9T) for some fixed e G (0,T/2). Then, the foUowing 
identity 



|z(x,0)| 2 cte 



d 
dt 



rj 2 (t)\z{x,t)\ 2 dx dt 



2K 



T Jf2 



7] 2 (t)z(x,t)z'(x,t)dxdt +2 / / 7/(t)r/(t)|z(x,i)| 2 (ixcft 



holds true for every z G H l (—T, T; L 2 (Q,)). Applying Young's inequality, this entails 

\z(x,0)\ 2 dx < (s + 2HVH00) / \z(x,t)\ 2 dxdt + s' 1 \z (x,t)\ 2 dxdt, 

Jq Jq 

for each s > 0, so the result follows by taking = 2||r7 / || 00 . 
Having said that, we are now in position to prove the: 

Lemma 3.3 There exist two constants C4 > and S4 > such that the following estimates 

2 



s 2 Cas / e 2 ^ 



n 



H 2 + |A| 2 ) dx - e 2s ^ (|VH 2 + |VA| 2 ) dx^j < f^3 k . 



lis) 



and 



hold true for k = 1,2, and s > S4. 

Proof . By applying Lemma I3T21 for z = e^U^i, we get that 

,,2 



2svo |U M (x,0)| 2 cix < s / e 2sip (s 2 \U kA {x,t) 2 \ + \\J kt2 {x,t)\ 2 ^dxdt, 



Cs l 



provided s is large enough. In light of (13.6I )- (13.7I ). this entails 



Cs 2 / e 2sw 



□ 



(3.10) 



C-,3 I e 2s ^ (|V//| 2 + |VA| 2 ) dx-J2 I e 2sLpo (\d a ^\ 2 + |d a A| 2 ) dx < ]T 3 fcJ (s), (3.11) 

JCl V ' H=2^ ' j,k=l 



curl(^B^) da;<3fc,i(s). 



Further, taking into account that ixBq G Hq (Q) 3 since /x vanishes in u and div Bq = 0, we have 



(3.12) 



/xBg dx-s 2 I e 2s ^° \V fi\ z dx < s 2 I e 2s ^° curl (^Bg 
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by (13.91) . whence 



Cs 3 [ e 2s ^° fiB^dx-s 2 [ e 2sw \Vfi\ 2 dx <3k,i( s ), 
Jn Jn 



(3.13) 



from (13.121) . Similarly, by arguing as above with z = e sip Vk,i instead of e^Ufc,!, we find some constant 
C > for which Cs 3 J Q e 2sipo j AD§ | 2 dx - s 2 j n e 2sipo |VA| 2 dx can be made smaller than the right 
hand side of (l3~T3l) by taking s sufficiently large. This, CD3> and (l3~T3l) entails (l3~T0l ). 

We turn now to showing (13- lib - To do that we apply Lemma [3721 with z = e^SjUj^i, i = 1,2, 3, 
getting 

C y e 2sw |a i [/ fc ,i(x,0)| 2 dx<s y e 2 ^(|VU fe ,i(x,t)| 2 + s- 2 |VU fci2 (x,t)| 2 )dx^, 
for s large enough. This yields 

C [ e 2s ^ d iC url(/iBg) ' dx < V3ij(s), 



(3.14) 



i=i 



by (I3.6I) - (I3.7I ). Further, bearing in mind that div Bq = and using that (di^B^ € Hq ($1) 3 , we obtain 



Cs [ e 2s v° (di/j)B% 2 dx- V / e 2sw |9 Q //| 2 dx < / e 2s ^° curl ((^)B§) 
Jn \ a \=2 



dx, 



(3.15) 



by (13791 . Moreover, as curl ((difi)B^) = djcurl (>B§) - ^curl ($.B§) -V/tx and /i G flj(n), 
we have by applying the Poincare inequality 

C [ e 2sipQ curl ((a^)Bg) 2 dx < [ e 2sLp0 ^curl (//B§) 2 + / e 2 ^°|V^| 2 dx, (3.16) 

in Jn 



hence 



(a^)B^ 2 dx- V I e 2s ^°\d a n\ 2 dx< I e 2sLp0 diCurl(juBg) 

Jn 



dx, (3.17) 



l<|a|<2 



by substituting the right hand side of (I37T61) for J* n e 2sw |curl ((^/i)B§)| 2 dx in (l3~T5l) . Putting (l3~T4l 
and (13.171 ) together, and summing up the obtained inequality over i = 1, 2, 3, we end up getting that 



Cs I e 2sv3 °\Vfi\ 2 dx- ^ / e 2 ^°|9V| 2 dx< 3*j(«). 

1<H<2^ 3,fe=l 



(3.18) 



Here we used again (13 - 8b - Finally, by arguing as before with z = e SLp diVk,i, i = 1,2,3, instead of 
^jUfc i, we get that (13.181 ) remains true with \i replaced by A. This completes the proof of (13- lib - □ 

Finally, we establish the: 

Lemma 3.4 There are two constants C5 > and S5 > such that we have 

2 

C 5 ! e 2s ^ (|9>| 2 + |«9 Q A| 2 ) dx - f e 2s ^ (|V/i| 2 + |VA| 2 ) dx < ]T 3k», 

\a\=2^° ^ n j,k=l 

for all s > S5. 
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Proof . In light of the two following basic identities \Jk,i(x, 0) = /^curl Bq + V/z x Bq and Vk,i(x, 0) 
— /j,cut\Dq — VA x Dq, k = 1, 2, arising from (13.6I ). we have 

\ 



/ Ui,i(x,0) \ 
-Vi,i(x,0) 
U 2j i(x,0) 
\ -V 2il (x,0) J 



, with „4(x) 



hence 



K{x) [ V *" ) = di A ( " 
\ V&A / \ A 



VA 



/ curlBj 
curl 
curl B 2 , 



curlD 2 , J 



( diV 1A (x,0) \ 
-8iV 1A (x,0) 
5iU 2 ,i(x,0) 



for every i = 1, 2, 3. From this and (11.111 ) then follows that 



£ (|a>| 2 + |a a A| 2 ) <c ^(|vu M (x,o)| 2 + |vv M (x,o)| 2 )+ ^ (>>| 2 + |d Q A| : 



|a|=2 



,fc=l 



|a|<l 



(3-19) 

Further, by multiplying (13.191 ) by e sipo , integrating over Q,, and upper bounding f n e 2 ^ ) VW^iXx, 0)| 2 <ix, 
with the aid of Lemma [3721 we find out that 



C [ e2SLp ° {\d a ^\ 2 + \d a \\ 2 )dx- I e 2s ^° (|V^| 2 + |VA| 2 ) dx 

|a|=2 

- s ( I e 2sip \VW ktl (x,t)\ 2 dxdt + s~ 2 [ e 2sLp \VVJ k2 {x,t)\ 2 dxdt) . (3.20) 
k=12 \JQ JQ J 

Here we took advantage of the fact that both fx and A belong to Hq(Q,) in order to get rid of the integral 
In e2sLpo (jH 2 + | A| 2 ^ dx by applying the Poincare inequality. Evidently the result now follows from 
(137201 and Lemma [2721 □ 



3.3 Completion of the proof of the main result 

In light of (13.111 ) and Lemma [3741 we may find C > such that 



r 2 
C J2 e 2s ^ (|9V| 2 + \d a X\ 2 ) dx<^ 3 fc)i ( s ), 



|a|<2""" 

upon taking s sufficiently large. Moreover, due to (I3.ll ). we have 



f e 2 ^ (|h feji | 2 + |Vh fc ,/) dxdt 

J Q 



< 



(B5,D£ 



C 3 (-T,T;W 2 -°°(n)) \ ^ 
|a|<2 



J e 2sv (\d a fi\ 2 + \d a x\Y dxdt^j , 



(3.21) 
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for every (J, k) G {1, 2} 2 , from where we get 

V ( C f e 2s ^° (\d a fi\ 2 + \d a X\ 2 ) dx- I e 2sip ( \d a fi\ 2 + \d a X\ 2 ) dxdt] 

2 

by combining (13.71) and (13-2 lb . Further, by recalling (12.71 ) and (12.101) we see for each a E {1, 2, 3} 2 with 
\a\ = 2, that 

! e 2sv (\d a n? + \d a \\ 2 ) dxdt = j e 2s ^° (|<9 Q /x| 2 + |9 a A| 2 ) G s (x)dx, (3.23) 

where 

/•T /-T 



:= / e -2»(w)W-ip(*,0)^ < / e- 2s{1 - aW) dt := g(s), a(t) := e" 7/3 * . (3.24) 

J-T J-T 

As lim s _ 5 . +00 g(s) = by Lebesgue's Theorem, we thus obtain from (13.22l )- (13.24l ) that 

r 2 

CJ2 ^ (l^l* + \d a X\ 2 ) dx<Y, (^-^(Wfcj) + s 3 e M ° s ||W fcJ |& 1(Q) ) , (3.25) 

|a|<2^ n * j,k=l 

for s sufficiently large. 

Furthermore, bearing in mind that xq E M 3 \ O, we notice from (12.71) and (12.101) that 

<Po(x) > mine 7 ^-^! 2 ^ ) > d t > d , x E H, (3.26) 
where do is defined in (12- 1 4b - From this and (13.251 ) then follows that 

r 2 

(\d a fx\ 2 + |d a A| 2 ) dx < e Cs ^(Wjfej) + s 3 e- 2( - dl - do >M, (3.27) 

[a[<2^ Q ' j,k=l 

where M is the same as in (11.121) . and 

#(W) = J (|V T V r | 2 + \V T XJ U \ 2 + |U;| 2 + |V;| 2 + (jU^I 2 + |V T | 2 )) dadt, W = (U,V). (3.28) 
In view of (f3T26b . (|3T27T )- (|3T28T) then yields 

Y I (\d a n\ 2 + \d a \\ 2 )dx<c{ ]T^(w feJ )) , 

\a\<2 Jfl ' \j,k=l J 

for some k G (0, 1), proving Theorem [T] 
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